Oscillatory angular dependence of the magnetoresistance in a topological insulator 
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The angular-dependent magnetoresistance and the Shubnikov-de Haas osciUations are studied in 
a topological insulator Bio.91Sbo.09, where the two-dimensional (2D) surface states coexist with a 
three-dimensional (3D) bulk Fermi surface (FS). Two distinct types of oscillatory phenomena are 
discovered in the angular-dependence: The one observed at lower fields is shown to originate from 
the surface state, which resides on the (211) plane, giving a new way to distinguish the 2D surface 
state from the 3D FS. The other one, which becomes prominent at higher fields, probably comes 
from the (111) plane and is obviously of unknown origin, pointing to new physics in transport 
properties of topological insulators. 

PACS numbers: 71.18.+y, 72.25.-b, 73.25.-|-i, 85.75.-d 
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The three-dimensional (3D) topological insulator hosts 
a novel surface state whose existence is guaranteed by a 
topological principle characterized by the Z2 invariant [l|- 
1^ , which makes it an interesting playground for exploring 
the role of topology in condensed matter l3-[20j . The key 
feature of the topological insulator is that an odd num- 
ber of massless, spin-helical Dirac cones comprise the 2D 
surface state, whose gapless nature is protected by the 
time-reversal symmetry. Such an intrinsically conduct- 
ing surface is supported by an insulating bulk, where 
the energy gap is created by a strong spin-orbit cou- 
pling. To understand the macroscopic properties of this 
novel surface state and to investigate the possibility of 
its device applications, transport studies are obviously 
important. However, in actual samples of topological in- 
sulators, there is always some bulk conductivity due to 
residual carriers, and separating the contributions from 
2D and 3D states turns out to be challenging [l3l4l5| |. 

Recently, we observed fi3\ strong de Haas- van Alphen 
(dHvA) oscillations in high-quality bulk single crystals of 
Bii^ajSb^: alloy in the "insulating" regime (0.07 < x < 
0.22), which is the first material to be known as a 3D 
topological insulator 0, Q. The dHvA oscillations sig- 
nified a previously-unknown Fermi surface (FS) with a 
clear 2D character that coexists with a 3D bulk FS [isj. 
Since Bii-jjSbj: is a 3D material, the observed 2D FS is 
naturally assigned to the surface. In the present work, 
to specifically probe low-dimensional properties of this 
material, we have extended our study to the angular de- 
pendence of the magnetoresistance. This method was 
successfully applied to the studies of quasi-2D organic 
conductors in the 1980s, resulting in the discovery of 
the celebrated angular-dependent magnetoresistance os- 
cillations (AMRO) |2l| - (23{ , and later extended to quasi- 
one-dimcnsional conductors [iJ,!!^. In this Letter, we 
present a detailed study of the magnetotransport in a 
Bio.9iSbo.o9 single crystal, where we have found oscilla- 
tory angular dependences in both the resistivity pxx and 
the Hall resistivity pyx- The oscillatory magnetoresis- 
tance (MR) is obviously different from known AMRO 



[25|, and it is comprised of two components: one is a pe- 
culiar manifestation of the 2D FS, and the other appears 
to be a fundamentally new effect. The former provides a 
new way to distinguish the contributions of the 2D FS, 
while the latter points to new physics in topological in- 
sulators. 

High-quality Bii_2:Sb2, crystals were grown from a sto- 
ichiometric mixture of 99.9999% purity Bi and Sb ele- 
ments by a zone melting method [13[ . The resistivity was 
measured by a standard four-probe method on a rectan- 
gular sample. In this Letter, we focus on a a; = 0.09 sam- 
ple in which the current was directed along the Ci axis. 
Continuous rotations of the sample in constant magnetic 
fields was used to measure the angular dependence of 
the MR within two main crystallographic planes (C3-C2 
and C3-C1). To observe the SdH oscillations, Pxx{B) and 
Pyx{B) were also measured by sweeping B between -fl4 
and -14 T along a set of magnetic-field directions. 

Let us start by presenting our main observations. Fig- 
ure 1(a) shows the angular dependences of the trans- 
verse MR measured in magnetic fields rotated within the 
trigonal-binary (C3-C2) plane perpendicular to the cur- 
rent. The magnetic-field strength, ranging from 0.01 T 
to 16 T, were kept constant during each rotation. As 
clearly seen in Figs. 1-3, pronounced oscillations of both 
Pxx and Pyx as a function of the rotation angle 9 appear 
for _B > 1 T. Two different types of oscillations can be 
distinguished: The first type consists of oscillations ap- 
pearing at lower fields (Figs. 1 and 2), while the second 
one becomes prominent at higher fields (Figs. 1 and 3). 
These data were taken immediately after the sample sur- 
face was refreshed by a chemical etching [2^. 

Figure 2 magnifies the "low- field" oscillations of Pxx{0), 
which become prominent above 1 T. These oscillations 
are symmetric with respect to the C3 axis {9 = 0°) and 
show a clear tendency of shifting the peak positions closer 
to the center with increasing magnetic field. It is impor- 
tant to notice that this strong field dependence is very 
different from ordinary AMRO [2^ where the angular 
positions of peaks are determined by commensurate mo- 
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FIG. 1: (Color online) (a) Angular dependences of Pxx 
measured in the trigonal-binary (C3-C2) plane in constant 
magnetic fields. The dashed line is the MR background 
~ (f — cos^ 6) for _B = 14 T. Inset depicts the measurement 
geometry, where 6 = Q° corresponds to the C3 axis, (b) Angu- 
lar dependences of pyx measured in the same conditions as in 
(a). The dashed lines are the expected angular dependences 
of the Hall effect plotted for all magnetic fields (0.1, 0.25, 0.5, 
1, 2, 3, 4.5, 6, 8, 10, 11, 12, 13, 14, 15, and 16 T). 



tions of electrons for certain field directions on warped 
FSs [2J, |25[ and, thus, are field- independent. Peculiarly, 
the amphtude of observed oscillations decrease with in- 
creasing field above ~8 T and disappear at higher fields. 
It turns out that the angular positions of the outer-most 
peaks (with respect to 6 = 0°), which are marked by ar- 
rows in Fig. 2, follow a simple geometrical law (inset of 
Fig. 2) related to the response of the 2D electrons resid- 
ing on the (211) plane perpendicular to the Ci axis, as 
will be discussed later. Note that, because of the crys- 
tal symmetry, there are three equivalent Ci axes in the 
plane perpendicular to the trigonal (C3) axis; hence, two 
of the three equivalent 2D FSs are at ±30° to the rotation 
plane. The pxx{0) measurements of the same sample in 
magnetic fields rotated within the trigonal-bisectrix (C3- 
Ci) plane (not shown) reveal essentially similar angular 
dependences. 

In contrast to these "low-field' oscillations, the "high- 
field" angular dependence of the transverse MR is devel- 
oping on a smooth field-dependent background coming 
from the anisotropy of Pxx{B) along the different axes. 
An example of its fitting for S ~ 14 T is shown by the 
dashed line in Fig. 1(a). Because of a large MR back- 
ground in strong magnetic fields, only largest peaks in 




FIG. 2: (Color online) "Low-field" angular dependences of 
pxx shown in Fig 1(a). Curves are shifted for clarity. Arrows 
mark the positions of the first peaks, which appear as the field 
direction is rotated from the C2 {C2) to C3 axis. Inset shows 
the dependence of 61 on B plotted as sin Si vs which 
gives evidence for the 2D-FS origin of the peaks. 



Pxx{0) can be clearly seen in the raw data [Fig. 1(a)]. 
Subtracting the background would help to reveal smaller 
peaks, but it is not necessary, because the situation is 
much simpler in the Hall data as shown below. 

Figure 1(b) shows the Pyx{0) data, which also show 
pronounced angular-dependent oscillations at high fields. 
The "background " for PyxiP) is simply the angular de- 
pendence of the Hall effect, RhB cos 9, where Rh is the 
Hall coefficient. As can be clearly seen in Fig. 1(b), low- 
filed Pyx{0) data follow this expected angular dependence 
very closely (we use Rh = —37 cm^/C, obtained from the 
Hall measurements), and the large deviation from this 
simple behavior is observed only in magnetic fields above 
10 T. Figure 3 shows "pure" oscillations in Apyx{9) af- 
ter subtracting the RhB cos 9 contribution from Pyx{9). 
One can clearly see a set of peaks, which are marked by 
short vertical ticks in Fig. 3. They are symmetric with 
respect to the C3 axis and show a rather complicated 
magnetic-field dependence. 

To understand the origin of the observed angular- 
dependent oscillations, detailed knowledge of the FS is 
useful. We therefore measured the SdH oscillations on 
the present sample for a series of magnetic-field direc- 
tions in two high-symmetry planes, as shown in Fig. 4(a) 
for the C3-C1 plane. Those SdH oscillations consist of 
three fundamental frequencies, which we call /i, /2, and 
/3. How these frequencies vary in the C3-C1 plane is 
shown in Fig. 4(b), which is essentially consistent with 
what was observed with the dHvA oscillations [13] and 
indicates the presence of both 2D and 3D FSs. In the 
present sample, the radius of the circular 2D FS (given 
by /i) is kp — 4.15x10'^ cm^^ which is the same as in 
Ref. 13, while the 3D FS is somewhat smaller, the semi- 
axes of the three ellipsoidal 3D FSs (given by /2 and 




FIG. 3: (Color online) "High-field" angular-dependent oscilla- 
tions of Apya: obtained by subtracting the expected RhB cos 6 
behavior from the py^ data shown in Fig. 1(b). Ticks mark 
the positions of distinguishable peaks. 



FIG. 4: (Color online) (a) SdH oscillations measured in the 
C3-C1 plane. For better visibility, dpxx/dB is calculated after 
a smooth background MR is subtracted, and the curves are 
shifted, (b) ^-dependences of the three fundamental frequen- 
cies obtained from the Fourier transform of the SdH oscilla- 
tions. Inset shows a schematic picture of the bulk Brillouin 
zone, in which the 2D FS is crossed by a plane perpendicular 
to the magnetic-field direction. 



/s) are a = 1.2x10^ cm~^, b = 8.5x10^ cm~^, and c = 
4.2x10^ cni~^ Note that the 2D FS resides on the plane 
perpendicular to the Ci axis, while the 3D ellipsoids are 
located at the L points of the bulk Brillouin zone, with 
the longest axes in the binary plane and tilted from the 
Ci axis by ^5°. The divergence of /i at 6' = 0° seen in 
Fig. 4(b) is of geometric origin and is a hallmark of the 
2D FS, as can be seen in the inset of Fig. 4(b). 

The analysis of the SdH amplitude gives us the cy- 
clotron mass Trie and the scattering time t. For the field 
direction very close to the C3 axis, the SdH oscillations 
coming from the 2D and 3D FSs can be easily sepa- 
rated, and their fits to the standard Lifshitz-Kosevich 
theory ^ at 9 ^ 3° yield rric/me of 0.13 and 0.033 for 
the 2D and 3D FSs, respectively {nie is the free elec- 
tron mass). Note that the relatively large rric for the 
2D FS originates from its characteristic angular depen- 
dence rric — rrfici sin0, where mP^ = 0.0057me is for the 
orbital motion in perpendicular fields {B \\ Ci), meaning 
that the 2D electrons are actually extremely light. Cor- 
responding Dingle plots (not shown) give Dingle temper- 
atures Td of about 5.5 K and 7 K for the 2D and 3D 
FSs, which imply the mean free paths of l^^ w 150 nm 
and « 16 nm, respectively [28j . This difference is not 
surprising, if one remembers the Dirac nature [29I Isoj 
of the 2D surface states and that they are topologically 
protected against spin-conserving backscattering. 

We now discuss the origin of the "low-field" part of 
the angular-dependent oscillations in Bio.91Sbo.09- A key 
insight comes from the recognition that the 3D FSs in 
our sample enter the quantum limit (where all the elec- 
trons condense into the first Landau level) with magnetic 
fields of only 1-3 T for any direction, while the necessary 
magnetic field to bring the 2D FS into the quantum limit 
diverges as the field direction is rotated toward the C3 



axis. This means that, in sufficiently high magnetic field, 
the 3D electrons will always remain in the quantum limit, 
while the 2D electrons can always be brought from the 
quantum limit back to the regime where higher Landau 
levels are populated. This will lead to the oscillating be- 
havior of pxx due to the 2D FS alone when the magnetic 
field is rotated toward the C3 axis. Remembering that 
the resistivity oscillates as p^x ^ cos[27r(i^/i?off +7)] with 
Bcff = B sin 0, a maximum in pxx occurs when F/ Sgff +7 
= n with integer n. Hence, one can understand that the 
Fermi level is crossed by the n-th Landau level of the 2D 
FS at the specific angle On given by 



arcsm 



(1) 



where F is the frequency and 7 is the phase of the SdH os- 
cillations. According to this model, the outer-most peaks 
marked by arrows in Fig 2 are coming from the crossing 
of the Fermi level by the first Landau level (n = 1) at a 
given magnetic-field strength, so their angles should cor- 
respond to Oi. Actually, the inset of Fig. 2 demonstrates 
that Eq. (1) fits the data reasonably well. If F is as- 
sumed to be 0.65 T based on the SdH result [3l|, the 
fit in Fig 2 yields 7 of 0.25, suggesting a non-zero Berry 
phase associated with the 2D FS 32l . 



Finally, let us discuss the origin of the "high-field" 
peaks shown in Figs. 1 and 3. Here, an important fea- 
ture is that the amplitude of these peaks weakens as the 
magnetic field is rotated away from the C3 axis, which 
is somewhat reminiscent of the behavior of the ordinary 
AMRO in quasi-2D systems if the co nduction planes lie 
perpendicular to the C3 axis 



21|-|23|. 



Note that this is 
different from what is observed for the "low-field" oscilla- 
tions, which are originating from the surface state resid- 
ing on the (2li) plane (_L Ci). Thus, it is probable that 
the "high-field" oscillations are coming from the (111) 
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FIG. 5: (Color online) (a) Angular-dependent oscillations of 
pyx within the C3-C2 plane in the constant field _B = 14 T 
measured at 1.5, 7, 15, 25, and 40 K. Inset shows oscilla- 
tions of iS.pyx after subtraction of the RhB cos 9 term, (b) 
Angular-dependent oscillations of p^x within the C3-C2 plane 
measured in the same conditions as in (a). 



plane (_L C£i, where surface states are seen in photoe- 
mission [THlOj and tunneling 11 1 experiments. Another 
distinguishable feature of the "high-field" oscillations is 
that they survive up to rather high temperatures as can 
be seen in Fig. 5, where temperature dependences of 
both Pyx{0) and Pxx{d) are shown. Even at 40 K there 
are still visible traces of oscillations (see the inset of Fig. 
5), while the "low-field" SdH and angular-dependent os- 
cillations disappear around 20 K. In spite of some simi- 
larities to the quasi-2D AMRO, the peak positions of the 
"high-field" oscillations apparently shift with the mag- 
netic field, which is not expected for the ordinary AMRO. 
Moreover, the existence of a finite coupling between con- 
duction planes is essential for the quasi-2D AMRO [25| . 
but there is no such inter-plane coupling for the surface 
states as long as the crystal is thick enough. Therefore, 
the observed "high-field" angular oscillations are a new 
phenomenon apparently specific to topological insulators, 
and it is possible that they are associated with a coupling 
between the surface and the bulk states. 

In conclusion, we found two different types of oscilla- 
tory phenomena in the angular-dependent MR of a topo- 
logical insulator Bio.91Sbo.09. One type is observed at 
lower fields and provides a new way to distinguish the 2D 
FS. The other one, which becomes prominent at higher 
fields, probably comes from the (111) plane and points 
to new physics in transport properties of topological in- 
sulators. 
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